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We discuss masses of stellar black holes found in binary systems and errors in their deter- 
mination. The observed mass distribution has a broad shape within the range 4 — IQMq 
without visible concentration to some preferred value. On the other hand, the initial black 
hole mass funciton as inferred from observations of luminous X-ray sources in other galaxies 
shows a power-law form. We show that both dynamically obtained black hole mass function 
and derived from X-ray observations can be made consistent in the frame of the hypothesis 
of enhanced black hole evaporation on a RSII brane for reasonable values of the warp-factor 
of the Anti-de-Sitter bulk space. 



1 Introduction 



There is a growing number of dynamical mass determinations of invisible compact objects in 
binary systems that are believed to be black holes (BH) primarily for their large mass (> SMq) 
and lack of evidence for having solid surface (see e.g. Orosz^-, Cherepashchukl^ for recent 
reviews). The mass distribution of compact objects in binary systems shows that masses of 
neutron stars (NS) are concentrated within a narrow range M^s = (1 ~ 2)M0 while all stellar 
BH masses fall within a wide interval ~ 4 — IQMq with a "mass gap" 2-4 Mq (e.g. Bailyn et 
al. 1^,1^ and references therein). Indirect information on BH masses in binary systems can also 
be obtained from examining properties of X-ray luminosity function of the most luminous (with 
X-ray luminosity > 10^^ erg/s) X-ray point sources observed in external galaxies by X-ray 
satellites assuming these sources to contain stellar mass BH and their luminosities being at the 
Eddington level. While dynamical BH masses appear to have a nearly flat distribution within the 
4-16 Mq range dN/dM ~ M^ , the X-ray luminosity function evidences for a power-law BH mass 
distribution dN/dM ~ M'^-^--^ ^. Unless some selection effects shape the form of BH mass 
distribution in both cases, BH mass functions derived from dynamical and X-ray observations 



are very different. Here we liypotliesize tliat the above properties of the observed BH mass 
function can be explained assuming an initial power-law BH mass function [dN / dM)^ ~ M~'^ 
with a ~ —2... — 3 and secular mass decrease of BH mass which is possible in the framework of 
some modern theories of multidimensional gravity. 



2 Determination of black hole masses and their errors 



The detailed analysis of systematic errors in determination of BH masses from observations 
which can affect the observed BH mass distribution is given elsewhere 1^1 and we only briefly 
describe them here. 



2. 1 Dynamical BH mass determination in binary systems 

The basic information on the mass of the invisible component in a binary system is obtained 
from the mass function of the optical companion 

f^im) = ."^y;, = 1.038 X 10-^K3p(i _ ^2)3/2 

where is the semi-amplitude of the optical star radial velocity curve (in km/s), niy and 
are masses of the optical and invisible star (in solar units), respectively, P is the orbital period 
of the binary system (in days), e is the orbital eccentricity and i is the binary inclination angle 
with respect to the line of sight. The mass of the invisible companion as derived from the optical 
mass function is ^ 

m^ = Mm) + (2) 
V qj sm-^i 

where q = rux/rno is the mass ratio. The only value of fv{fn) immediately provides the lower 
limit rrix > fv{fn). More precise value of rux requires the knowledge of q and sini. 

a) Optical star form effects. 

If the mass ratio q<l (high-mass X-ray binaries, HMXB, e.g. Cyg X-1, LMC X-1, SS 433)), 
the system's barycenter lies within the optical star body and tidal distortions of the optical star 
form mostly influence the spectral line profiles which are used to determine the radial velocity 
curvel^. Another distortion of radial velocity curve and its semi-amplitude in HMXB is due 
to variable selective absorption of light by strong stellar wind of the massive 0-B optical star^^. 

In low-mass X-ray binaries (LMXB) g > 1 the system's barycenter is outside the optical 
star so the effects of the star form can be neglected. Stellar wind effects from low-massive A-M 
dwarfs in these systems are also insignificant. 

b) Stellar rotation effects. 

The mass ratio q is usually derived from rotational broadening of absorption lines in the 
optical star spectrum. Indeed, assuming synchronous axial and orbital rotation, one can 

/ l\2/3 

vs\ni = QAmK^q-^/^ [l + -j . (3) 

The determination of vsvai from usual analysis of absorption line profiles can not be made 
better than to within the 10-20% errors due to X-ray heating of the optical star atmosphere 
However, parameters q, i can be derived from orbital variabilit y of the absorption line profiles 
using new method proposed by Antokhina and Cherepashchuk^^ and Shahbaz^^. This allows 
the strong reducing of errors in the parameter q determination. 

c) Binary inclination angle effects. 

Uncertainties in the binary inclination angle i provide the largest error in the ma ss d eter- 
mination. The usual way of measuring i is from the ellipticity effect of the optical star^^. The 



main uncertainty here comes from model-dependent contributions of other emitting structures 
(gaseous stream, accretion disk) into the total optical or infrared variability of the system and 
is especially important (> 50%) in LMXB (in which most BH candidates have been discovered). 
The new method of determinat ion o f parameters g, i based on the analysis of the orbital variation 
of the absorption line profiles independent of contribution from other gas structures in 

the binary system. However, it can be applied at a very high spectral resolution of 50000-100000 
and can be realized only on very large telescopes. 

So, with the 90% certainty we can state that masses of stellar BH measured dynamically 
in binary systems span a wide rage 4-16 Mq without visible concentration to some value, i.e. 
dN/dM ~ M° with a ~ 0. 

2.2 X-ray luminosity function of ultra-luminous X-ray sources 

Ultra-luminous X-ray sources (ULX) are point-like sources with persistent X-ray luminosity 
Lx > 2 — 410^^ erg/s corresponding to a 1.4 Mq neutron star. They have been discovered by 
several X-ray satellites (mostly by Chandra., see e.g. '•^ and references therein) in other galaxies. 
Although some of them proved to be another galaxies seen through the observed galaxy disk (see 
^Ifor a recent identification), they constitute the most luminous population of X-ray sources in 
galaxies and based o n th eir X-ray spectral fitting can be identified with accreting stellar-mass 
BH in binary systems^^. 

A general a naly sis of X-ray source populations in nearby galaxies was recently carried out 
by Grimm et al. It revealed the universal power-law shape of the average X-ray luminosity 
function dN/dL^ ~ L~'^'^ in a wide range of X-ray luminosities 10'^^ — 10^^ erg/s with a steeper 
decline oc L~'^''^ for the most luminous sources > 2 x 10^^ erg/s. The universal shape of the 
X-ray luminosity function below (2 — 4) x 10^^ erg/s can be explained on the very general grounds 



by accretion in binary systems'^. The cut-off in the luminosity function at Lx > W^^ erg/s (i.e. 
for ultra-luminous X-ray sources) appears to be a general feature of X-ray source populations 
in galaxies (see e.g. a recent analysis of Chandra observations of NGC 5194 in^^ reflecting the 
presence of binaries with BH accreting at the Eddington limit. This allows us to assume that 
for BH in these systems dN/dM = {dN/dLx){dLx/dM) ~ ^ 

Uncertainties in this method of determining BH mass distribution are due to possible beam- 
ing of X-ray emission from binary systems (the value of the X-ray luminosity is usually obtained 
from the observed X-ray fluxes assuming spherical symmetry). At present, this uncertainty 
is hard to estimate, but similar shapes of X-ray luminosity functions of point-like sources in 
different galaxies and the very possibility of obtaining a universal power-law form for dN/dLx 
which is simply scaled by the star formation rate from galaxy to galaxy (as argued in^^ seem 
encouraging. 

3 Black hole mass function and enhanced black hole evaporation 

Clearly, the decreasing power-law shape of BH mass distribution is drastically different from the 
nearly flat form derived from dynamical BH mass determinations. First, we should note that 
if ULX are actually BH accreting at the Eddington limit, the mass accretion rates M in these 
binaries must be above ~ lO~^M0/year, i.e. such systems must be HMXB with ages < 10^ 
years. In contrast, most BH with dynamically determined masses reside in transient LMXB 
(X-ray novae^, 121 which are old systems with ages > 10^ years. 

Let us consider the distribution function f{M) = dN/dM of a population of sources with 
changing mass. Let the initial mass distribution be fo{M) within the mass range [Mmin, Mmax] 
and the law of the mass change be M(t). In the stationary case the evolution of mass distribution 




is described by the one-dimensional kinetic equation 

d 



f{M)M =MM) (4) 



dM 

so the stationary distribution function for M > will be 

Sm UM')dM' 

f(^M) = )^ ^ M<M^acc (5) 

while for M > M^ax the stationary distribution is independent of the initial mass function and 
is determined only by the mass change M: 

/aY'""" MM')dM' const 

/(M) = -^^ ' — = iir' 

For example, assuming the initial power-law mass function fo{M) cx M~"' yields for > 1 
/(M) ~ [1 _ (M/M„i„)-°'+i]/Af ~ 1/M at M > a few M^i„, and for < 1 /o(M) ~ 
[(M/Mmin)^"'^^ — 1]/M. Since BH mass in accreting binaries can never increase faster than 
M oc M, the stationary distribution never (for any a^) decreases steeper than 1/M. So the ob- 
served steep BH mass cut-off (propto in ULX apparently evidences that stationarity 
arguments for these systems are inapplicable due to short time of Eddington-limited accretion 
in these systems (typically ~ 10^ years <C M/M) so that BH masses cannot change significantly. 
Then we stay with the possibility that the observed steep BH mass function reflects the initial 
BH mass distribution in these binaries. 

In LMXB (in which most BH candidates are found) or for single BH (only two of them 
are suspected by gravitational microlensing experiments'^^ the situation can be quite different. 
Small mean accretion rates M ~ 10~^MQ/yr typical for these systems are insignificant for BH 
mass growth. Applying the same arguments as above, we could expect the initial BH mass 
function to have a nearly flat shape in this case. However, there is another possibility. 

Let us consider the hypothesis of enhanced BH mass evaporation recently put forward by 
Tanaka^^ and Emparan et al.l^. These authors consider the RSH brane world models I^Hl in 
which our Universe is localized on a 4D-brane embedded in an Anti-de-Sitter bulk characterized 
by the warp factor (curvature length) C. Such models admit macroscopic values of C up to 0.1 
mm so as not to contradict the existing laboratory measurements^^. In this setup, the authors 
El El speculate that no stationary macroscopic BH can exist on the brane due to enhanced 
evaporation into (virtually unobservable) CFT-modes (KK-gravitons). This process has not 
yet been calculated properly so this hypothesis remains highly speculative (see e.g. Casadio 
2003!^ who obtained different results). Nevertheless, we try to see what happens should the 
macroscopic BH evaporation actually exist. 

In this hypothesis, the evaporation time r of a macroscopic BH localized on the brane (i.e. 
with the horizon size > C) is shorter than the classical Hawking evaporation time ~ tpi{M /mpif 
by a huge factor 0{C/lpi)'^ (subscript PI stands for the Planck units) 

M \3 /Ifmml^ ^ 



So the mass of an isolated BH decreases as M M~^. It is easy to show that for C ~ 
oof~^ — oo/~^ mm the BH evaporation rate exceeds the mean accretion rate in LMXB for 
M ~ IOMq. So the present bounds C < f.oo mm do not contradict this assumption. 

Now for old BH in low-mass X-ray binaries we can find the stationary shape of BH mass 
distribution function. For evaporation M < so the integration of Eq. Q) yields 

/(M) = i2=lAp^£, M>M_ (8) 
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Figure 1: Schematic view of the expected stationary BH mass function (the solid curve) with the initial power-law 
form (dN/dM)i ~ M~"» (the dashed line) obtained in the enhanced BH evaporation model on the RSII-brane. 
The mass Mo corresponds to a minimal BH mass that can evaporate in the Hubble time. 



and 

nM^.l^tZJfni^.^, M<M„.,. (9) 

Assuming as above /o(M) ~ M""' with < 1 we have /(M) ~ M""'+3 f^j. j^j > j^^.^ 
and f{M) ~ for M < Mmin (we always assume M <C M„iax)- This is illustrated in Fig. 1. 
Clearly, if ~ 2, as can be inferred from ULX observations, the stationary BH mass distribution 
could be made much flatter by BH evaporation. In addition, the steep decrease of the stationary 
BH mass function (/(M) ~ M^) below Mmin can be identified with apparent deficit of BH with 
masses below 4 Mq. 



4 Conclusions 



The steadily growing number of mass determinations of stellar BH allows us to construct and 
study BH mass function. It has apparently broad and near flat shape within the interval 4— 16Mq 
without visible concentration at some value, which is opposite to what is observed for NS masses 
{Mjsfs = 1 — 2Mq). Possible systematic errors in BH mass determination, which mostly rely 
on the radial velocity curve of the optical companion, can be reduced in future high-resolution 
spectroscopic observations. 

BH mass function can also be inferred from observations of luminous X-ray binary systems 
(with Lx > 10'^^ erg/s) in other galaxies. The current Chandra observations suggest a steep 
decrease {dN/dLx ~ L^"^) in the X-ray luminosity function of point-like objects in galaxies. If 
these sources are accreting BH at the Eddington limit, the initial BH mass function should have 
a similar slope fo{M) ~ M~^, which disagrees with BH distribution found in galactic binary 
systems. 

We show here that in the framework of the enhanced BH evaporation hypothesis^^ which 
can in principle take place in some multidimensional gravity models (of RSII type), the initial 
steep power-law BH mass function can be made flatter to match the observed BH distribution 
in galactic binaries. This hypothesis also predicts low-mass BH (with M less than a few solar 
masses) should be very rare. Discovery of such low-mass BH would strongly limit this hypothesis. 



Concluding, we note that it is unclear at present how close to reality such multidimensional 
gravity models are, so any test of their astrophysical predictions is very desirable. We hope that 
future precision measurements of BH masses (both single and in binaries) by traditional and 
new (e.g., gravitational waves from coalescing binaries with BH) astrophysical methods can be 
potentially very interesting in this respect. 
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